Quantized Non-Abelian Monopoles on S 3 
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A possible electric-magnetic duality suggests that the confinement of non-Abelian electric charges 
manifests itself as a perturbative quantum effect for the dual magnetic charges. Motivated by this 
possibility, we study vacuum fluctuations around a non-Abelian monopole-antimonopole pair treated 
as point objects with charges g = in/2 (n = f , 2, ...), and placed on the antipodes of a three sphere 
of radius R. We explicitly find all the fluctuation modes by linearizing and solving the Yang-Mills 
equations about this background field on a three sphere. We recover, generalize and extend earlier 
results, including those on the stability analysis of non-Abelian magnetic monopoles. We find that 
for g > 1 monopoles there is an unstable mode that tends to squeeze magnetic flux in the angular 
directions. We sum the vacuum energy contributions of the fluctuation modes for the g = 1/2 case 
and find oscillatory dependence on the cutoff scale. Subject to certain assumptions, we find that 
the contribution of the fluctuation modes to the quantum zero point energy behaves as —R~ 2 ^ 3 and 
hence decays more slowly than the classical —R^ 1 Coulomb potential for large R. However, this 
correction to the zero point energy does not agree with the linear growth expected if the monopoles 
are confined. 

PACS numbers: 14.80.Hv, 11.15.-q 



I. INTRODUCTION 

There are many facets to the investigation of magnetic 
monopoles. The most straight-forward view is that mag- 
netic monopoles are classical solutions in certain gauge 
theory models p], 0| , and the possible realization of these 
solutions leads to important ramifications for particle 
physics and cosmology. Another facet is in the spirit of 
the original proposal by Skyrme to understand the proton 
as a classical solution 0, [j] . Then magnetic monopoles 
are themselves some versions of particles, and as sug- 
gested by the sine-Gordon and massive Thirring model 
equivalence [H, Q , magnetic monopoles are dual to "or- 
dinary" electrically charged particles whereby, in some 
region of parameter space, it makes better sense to view 
the magnetic monopoles as the fundamental particles and 
vice versa. A tantalizing correspondence of this kind ex- 
ists in grand unified inspired models where there is a 
one-one correspondence between the magnetic charges of 
stable monopoles and the electric charges of known par- 
ticles such as quarks and leptons @, [!]■ Also, the dual- 
ity between particles and magnetic monopoles has been 
shown in certain supersymmetric theories Q, substanti- 
ating the early conjecture in Ref. (loj . 

If magnetic monopoles are to be viewed as magnetic 
versions of electrically charged particles, perhaps it is 
simpler to study certain properties of electric charges in 
the magnetic sector. To follow this line of thought, we 
know that magnetic monopoles, like ordinary particles, 
can carry non-Abelian charges. Now since particles car- 
rying non-Abelian charges are thought to be confined, 
perhaps non-Abelian magnetic monopoles are also con- 
fined, and it may be simpler to understand confinement 
by studying magnetic monopoles instead of electric par- 
ticles. This argument motivated us to study quantum 
effects in the background of non-Abelian monopoles. 

Confinement occurs due to non-perturbativc effects in 



the non-Abelian electric backgrounds of particles. The 
non-perturbative effects are supposedly due to condensa- 
tion of magnetically charged objects. However, duality 
relates the electric and magnetic sectors with inversely 
related charges. Hence, if we work in the dual sector, 
the electric sources get replaced by magnetic charges 
(monopoles) and, if there is confinement, it would be 
due to the condensation of electric charges. We empha- 
size that the exchange of (weak, electric) and (strong, 
magnetic) is simply a hypothesis. However, within 
this hypothesis, magnetic monopoles should get confined 
by perturbative quantum effects in the electric sector. 
The quantum corrections to the magnetic field of the 
monopolc should be calculable using a perturbative ex- 
pansion in the electric coupling. 

It is to be noted that the effect we are looking for 
is distinct from the confinement of magnetic monopoles 
due to spontaneous symmetry breaking by the condensa- 
tion of an order parameter. We do not have extra scalar 
fields that take on a vacuum expectation value to break 
the non-Abelian symmetry and hence do not have any 
topological strings that confine the monopoles. Instead 
we are hoping to see evidence for a quantum string that 
confines the monopoles. 

Specifically, we aim to find the one-loop quantum con- 
tribution to the energy of a non-Abelian monopole and 
an anti-monopole as a function of their separation (see 
Fig. [1]). In this, we have to assume a background classi- 
cal field configuration for the monopoles and this is con- 
structed by linearly superposing the spherically symmet- 
ric classical solutions. If the quantum corrected energy 
grows faster than the separation, it would indicate that 
the choice of spherically symmetric magnetic field config- 
urations for each monopole is not preferred since, if all 
the magnetic flux is confined to a tube, the energy only 
grows linearly with separation. 

A monopole and an antimonopolc in R 3 is a compli- 



FIG. 1: A classical monopole and an antimonopole in R 3 and 
their magnetic field lines. The system has little symmetry, 
making it difficult to evaluate quantum fluctuations in this 
background. 



cated background in which to study fluctuation modes 
which then have to be summed up to find the quantum 
corrected energy (see though (HI). Hence we have cho- 
sen to study a monopole and an antimonopole placed at 
the antipodal points of an S 3 , as shown in Fig. [21 The 
classical magnetic field of the monopoles is then spher- 
ically symmetric, and mitigates the difficulty of finding 
the fluctuation modes. For the time being, we will sim- 
ply focus on the analysis on S 3 , and hope that a mapping 
to R 3 can be made at a later stage, or perhaps a similar 
analysis will be found to be feasible directly in R 3 . 

Our first task is to construct non-Abelian monopoles. 
The simplest model with such monopoles involves the 
symmetry breaking 

SU(3) -> [SU(2) x U(l)]/Z 2 

and the fundamental monopoles have non-trivial SU(2) 
and U(l) charge. The symmetry breaking is achieved by 
giving a vacuum expectation value to a field $ trans- 
forming in the adjoint representation of SU(3>). The 
monopoles have a regular core and are known quite ex- 
plicitly since they are basically SU(2) 't Hooft-Polyakov 
monopoles embedded within SU(3) \12\ . 

Analysis of fluctuations around monopoles arising from 
the SU(3) model is very complicated since it involves 
both scalar and gauge fields. Also, the structure of the 
core will be important in determining the fluctuation 
modes. So we have chosen to study non-Abelian Dirac 
monopoles. In this case, described more explicitly in 
Sec. mi the monopole core is singular but the external 
magnetic field for g = 1/2 is identical to the SU(2) part 
of the fundamental SU(3) magnetic monopole. Ignor- 
ing the core structure does mean that we need to impose 
boundary conditions at the location of the monopole by 
hand, an issue that we discuss in Appendix [Bl 

In solving for the angular part of the fluctuation modes 
we are aided by the early work of Wu and Yang [l3[ 
who found monopole scalar harmonics, Brandt and Neri 
who studied the stability of non-Abelian monopoles, 
Olsen et al [l5| and Weinberg , who constructed 
monopole vector harmonics (also |14 ]. The radial prob- 
lem in R 3 has been set up in Ref. [l5[ though with the 



FIG. 2: A monopole and an antimonopole on S 3 and the 
classical, spherically symmetric field lines. (One dimension of 
the S 3 has been suppressed.) 



limited goal of finding a bound state solution. Instead 
our analysis requires us to evaluate all solutions of the 
angular and radial problems on S 3 , in the background 
of an monopolc-antimonopolc pair (and not just a single 
monopole). 

To evaluate the one-loop quantum correction to the 
energy of minimal charge monopoles with a = 1/2 on S 3 , 
we follow the procedures developed in [I?) (see [l9| 
for more detailed expositions). We find all the cigenfrc- 
quencies of the fluctuation modes and then sum them up 
to an ultraviolet cutoff. We do not face some of the sub- 
tleties in the quantization of solitons since we are in a 
compact space and so all the fluctuation modes are dis- 
crete. However, we treat the monopoles as point-like (i.e. 
non-Abelian Dirac monopoles) and have no control on 
the renormalization of its mass for example. But we are 
primarily interested in the dependence of the quantum 
corrections to the energy as a function of the monopolc- 
antimonopolc separation, and assuming that the mass of 
the Dirac monopole has no dependence on the size of the 
S , this is given by the dependence of the sum of zero 
point fluctuations on R, which we can indeed calculate. 

We begin by describing the background of non-Abelian 
monopoles on S 3 in Sec. [Til This is followed in Sec. lIIII bv 
a derivation of the linearized fluctuation equations. Here 
we see that the non-Abelian gauge fluctuations split up 
into 4 types of excitations depending on whether they 
are uncharged or charged and whether they are scalar 
or vector in nature. The uncharged scalar and vec- 
tor sectors do not feel the non-Abelian monopole back- 
ground but do play a role in the zero point energy, and 
we construct all of them on S 3 . Charged scalar fluc- 
tuation modes in the background of a single monopole 
in R 3 have been considered by Wu and Yang [l3[ lead- 
ing to the Wu-Yang spherical harmonics. Here we also 
solve for the radial dependence and find the full set of 
charged scalar harmonics on S 3 in the background of 
the monopole-antimonopolc. Similarly vector spherical 
harmonics in the background of a monopole in R 3 have 
been constructed in Refs. [3, HH, Here we find all 
the charged vector spherical harmonics on the monopole- 
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antimonopolc background in S 3 , including the radial de- 
pendence of the eigenfunctions. In Sec. HVI we have sum- 
marized the technical results for the derivation of the 
eigenvalues and eigenfunctions that are shown in detail 
in the Appendices. For the g = 1/2 monopole, the sum 
over vacuum fluctuations is done in Sec. IVIII We close 
the paper with discussions in Sec. I Villi 



II. MONOPOLE BACKGROUND 

We consider a Yang-Mills field on a three sphere S 3 of 
radius R. The S 3 metric in polar coordinates (x, 0, and 



ds 2 = -dt 2 + R 2 



d X 2 + sin 2 X 



d9 2 + sin 2 6d<p 2 I (1) 



The evaluation of the fluctuation modes is quite tech- 
nical and so it is helpful to summarize our findings qual- 
itatively at this stage for the reader who may not be 
interested in details. As far as the stability and possible 
confinement of non-Abelian monopoles is concerned, the 
charged vector modes are of particular interest. These 
vector modes interact with the monopole background via 
their charge and also by the Zeeman coupling of their spin 
to the background magnetic field. In addition, they expe- 
rience the curvature of the S 3 . If g > 1, where g = eg /Aw 
and g is the magnetic monopole charge, we find that there 
is an unstable mode in the spectrum of fluctuations, im- 
plying that the spherically symmetric background is un- 
stable. This is consistent with the result in [l4| but we 
also find the unstable mode explicitly and show that it 
tends to squeeze the monopole flux in the angular di- 
rections. For g = 1/2, however, this classically unstable 
mode is absent. In this case, we find the eigenvalues and 
eigenfunctions for all the fluctuation modes. We then 
evaluate the contribution of these quantum zero point 
fluctuations to the energy of the monopole background, 
as compared to the zero point energy of the trivial back- 
ground, by summing up the energies of all the individ- 
ual modes. If we sum up to T c modes, we find that 
this contribution to the vacuum energy oscillates as a 
function of r c . Averaged over the oscillations the con- 

1 /3 

tribution is proportional to — T c / R- If we impose an 
ultraviolet cutoff on the momentum, this corresponds to 
a — i?~ 2 / 3 dependence of the vacuum energy on the ra- 
dius of the S 3 which is directly related to the separation 
of the monopolc-antimonopole pair. 



Additional contributions to the vacuum energy will 
arise from renormalization of the parameters, in this case 
the mass of the monopoles. We have taken our monopoles 
to be point-like and hence do not have control over these 
extra contributions. In any case, the —R~ 2 / 3 dependence 
is to be contrasted with a +R growth in the energy if the 
monopole and antimonopolc arc confined by a flux tube. 
Hence the energy of the spherically symmetric contribu- 
tion does not become larger than the energy of a flux 
tube for some critical R and the spherically symmetric 
configuration is of lower energy. So we conclude that the 
g = 1/2 monopoles on S 3 are not confined, at least for 
reasons of vacuum energy. However, the quantum zero 
point energy may still turn out to be larger than the clas- 
sical —1/R Coulomb energy for monopolc-antimonopole 
separation larger than some critical value. 



The Yang-Mills equations of motion are 



M c 



(2) 



where F M „ = d^Au — d v A^ + ie [A^, A v }. The curvature 
of space is taken into account via the covariant derivative 
V M which acts on vector fields as V \A V = d^A" +T^ a A a . 
For simplicity we will consider only SU(2) Yang-Mills 
fields in this paper. The generalization to SU(N) is 
straightforward. 

An exact solution to the Yang-Mills equations (0) is 



1 



1 — cos 

AlT 



us} 



< 6 < it 



(3) 



where the direction in group space, and q is the 

magnetic charge. Note that the gauge potential is not 
defined at = n. As discussed in pjjj, we also need to 
define the gauge potential on a second coordinate patch 



1 + cos 

47T 



,(3) 



< < 7T 



(4) 



The gauge fields in the region of overlap of the two co- 
ordinate patches are related by a gauge transformation 
that is well-defined provided 



eg n 
9 ~ An~ ~ 2 



0,±1,±2,... 



(5) 



The background field strength can be calculated from 
the gauge potential and is 



sin0 



.(3) 



(6) 



and corresponds to the canonical monopole magnetic 
field 

x ?- 2 sm0 v An r 2 

where r = sinx- 

This exact solution corresponds to a monopole at the 
north pole (x = 0) of the S 3 and an antimonopole at 
the south pole (x = tt). The solution is singular at the 
location of the poles, as is the case for Dirac monopoles. 
This issue could be eliminated by considering regular 't 
Hooft-Polyakov monopoles P, [|J , but it would complicate 
the analysis significantly and would not change the long 
distance behavior of the non-Abelian system which is our 
primary interest here. 
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III. LINEARIZED YANG MILLS THEORY 

We now consider a linear perturbation a M around the 
background monopole solution A^. Perturbing Eq. ([2]) 
and keeping terms that arc first order in a gives 

gT { (fljAoa - r^Aaa) + r^ Q (9a a v - d v a x ) 
-Y x ( [a x ,A a ] + [A x ,a a ] 



riJ [a v ,A x ] + [A u ,ax] ) +2 [a M ,.F„ a ] 



+ (ie) 2 ( [A M , H„, oj] - [A Q , a„]] ) } 
+ (d a g^) ( d M a„ - r*„a A + «e a„] ) 
-.^(^r^K = . (8) 
We work in the fixed background gauge (2l| . 

<? Mi/ (5 M a„ - r^o A + *e a,]) = . (9) 

The equations of motion can be cast into a more friendly 
form by decomposing the perturbation in terms of its 
group components, as well as separating the time and 
spatial parts of the 4-vector, 



fc=±,3 

a« = (a« a (k) 



(fc) 



(10) 



Here k = ±,3 is the group index, <jw are the Pauli 
matrices with a^' = &^ ± i<5-( 2 ). With this decom- 
position, the equations decouple into 4 sectors: the un- 
charged time component which behaves as a scalar a t , 

the scalar charged component a{ , the uncharged spatial 
vector aS 3 \ and the charged vector a^ + '. The equations 
of motion can then be schematically written the following 
way: 



(-R 2 d 2 + A< 3 >) <4 3) = (11) 

R 2 d 2 + A( 3 >) a< 3 ) = (12) 

-i? 2 a t 2 + a«) 4 +) = o (13) 

-R 2 d 2 + A' + >) a (+) = . (14) 



Here A and A are the scalar and vector Laplacians re- 
spectively, and in the (+) direction the derivatives are 
gauged, di — > d% — ieAi. The vector Laplacians take 
into account the coupling to the curvature of the S 3 
sphere, via the Ricci scalar, and in the (+) direction 
there is also a Zeeman-likc coupling between the back- 
ground monopole and the perturbations. We have not 



written the equations for the (— ) scalar and vector com- 
ponents as they are simply the complex conjugates of the 
Eqs. (JT3J) and (JTU of the (+) direction. 

The gauge condition © decouples into the neutral and 
the charged sectors, but it does mix between the time 
component and the spatial vector. It can be schemati- 
cally written as 



-i?d t c4 +) +V(+>a(+) 



(15) 
(16) 



The explicit equations can be found in Appendix [X] The 
boundary conditions to which they are subject are dis- 
cussed in Appendix iBl 



IV. ENUMERATION OF EIGENMODES 

Constructing the solutions to the linearized Yang-Mills 
equations is a long yet subtle exercise in classical math- 
ematical analysis. The details are given in Appendix IB] 
here we simply enumerate the solutions and their associ- 
ated frequencies. At first we will simply list all the solu- 
tions to the equations for the fluctuations. In Sec. IIV Gl 
we will impose the gauge conditions and also eliminate 
modes that are pure gauge. Finally we will be left with 
the physical modes and these are summarized in Sec. [V] 



A. Neutral scalar sector 

In this sector we compute the eigenmodes of the scalar 
Laplacian on S 3 . The modes are derived in Appendix 
IBl by separation of variables and again in Appendix [C] 
by the use of group theory; the latter derivation exploits 
the four-dimensional rotational symmetry of the problem 
and shows how the solutions organize into multiplets that 
are irreducible representations of the so(4) algebra. The 
resulting solutions are: 



af ] =Yj M (e,d>)^j ! j(xV 



Here 



€n-J,j(x) = sin J (x)G; i Ij(cosx) 



(17) 



(18) 



where Yjm(&,4>) are the familiar spherical harmonics 
and G are the Gegenbauer polynomials (20| . The quan- 
tum numbers span the range n — 1,2,3,... and J = 
0, 1, 2, . . . , n and M = — J, . . . , +J. The energy of this 
solution is given by 



A = uR = \fn(n + 2) 



B. Charged scalar sector 



(19) 



In this sector we must compute the eigenmodes of the 
scalar Laplacian on S 3 that has been minimally coupled 



5 



to the gauge field of the monopole antimonopole pair. 
The solutions are obtained by separation of variables: 



4 +) =Wj M (0A)Cn-J, a (x)e i 



where 



ics [3 i 




(20) 



(21) 



are the Wu-Yang monopole spherical harmon- 
and £ n -j, a is given by 



£n-J, a = sin Q (x)G"^j(cosx)- 



(22) 



The quantum numbers satisfy n = g, g + 1, g + 2, ... and 
J = 9i9 + 1) 9 + 2, • • • n and M = — J ... J. The energy 
of this solution is given by 



A = yj(n - J + a){n - J + a + 2) 



(23) 



Note that a = J for g = and the charged solutions 
reduce to the neutral ones as expected. 



C. Neutral vector sector 

In essence we are looking for the eigenmodes of the 
vector wave equation on S 3 here. Since this is a vector 
equation we expect three sets of solutions corresponding 
to different states of polarization. One set of solutions, 
Eq. ([24)) , may be constructed simply by taking the gradi- 
ent of the neutral scalar eigenmodes, Eq. (fT7j) . The third 
set, Eq. (|26|) . is obtained by taking the curl of the second. 
The second set of solutions, Eq. (|2"5j) below, are obtained 
by looking for transverse solutions that have no radial \ 
component, as discussed in Appendix [Bl 
First solution: 

4 3) (i) = Y JM (e^)d x ^.jAx)e wt 

4 3) (1) - d e Y JM {e,cj>)^ n _j,j{x)e iwt 

af{\) = d^YjMie^^j.Axy^ 

A = y/n{n + 2) 
Second solution: 



(24) 



af\2) = -i^[^r /M (e,0)]sin X Cr l -J,j(x)e 1 ' 



sin 9 



af\2) = - sin 9[d e Y JM {0 A)} sin x£n-J,j(x)e iwt 

A = n+ 1 (25) 
Third solution: 

, „ 1 



a^(3) = J(J + 1K 



JM \ 



smx 



<n-jAxV' 



4 3) (3) = dgYjMie^^lsinxdn-jjix)}^ 1 
af(3) = d^Y JM {6,4>)d x [sin X 6w,j(x)] e 4 ' 



The quantum numbers span the range n = 1,2,3,..., 
J = 1, 2, . . . , n and M = — J, . . . , J for the second and 
third solution. For the first solution the allowed range is 
n = 1,2,3,..., J = 0, 1, 2, . . . , n and M = — J, . . . , J; in 
other words the value J = is also allowed. 



D. Charged vector sector 

There is a close parallel between the solutions in the 
neutral and charged vector sectors. With one excep- 
tion, the charged solutions can be obtained from the neu- 
tral solutions by the replacement of spherical harmonics 
Yjm with Wu-Yang harmonics Wjm, of J(J + 1) with 
J (J + 1) — g 2 , and the minimal coupling substitution 
d(j, — > drf, — iAfj, = where is the background vector 
potential. The exception is that the second solution in- 
volves an additional term proportional to g as explained 
in Appendix [Bl 
First Solution: 

4 +) (i) = WjMie^^d^n.jM^ 
4 +) (l) = d Wj M (e,<t>)Z n -J,a(x)e iut 

4 +) (i) = D4,WjM(e,<j>)Z n -j, a (x) eiu}t 



A = yj{n- J + a)(n - J + a + 2) (27) 
Second Solution: 
4 +) (2) = 

4 +) (2) = S m X ti(lJ(J+l)-g 2 }-^D^W + gid e W 



sin ( 



a< +) (2) = smxZ{-{J(J+l)-g Z }ism6deW + gD^W) 



A 



1-J + a 



(28) 



A 



71 + 1 



(26) 



For brevity we have written £„_j a (x) as £ and 
Wj M {6,(t>) asW. 
Third solution: 

4+>(3) = [J(J + l)-g 2 ]Wj M ^—Zn-j,ae iut 
A smx 

4 +) (3) = d e W j M {e,4>)d x [sin x^n-jAx)}e zwt 
a^ +) (3) = D^Wj M (0^)d x [sin X ^jAx)}e luJt 

A = n + l-J + a (29) 

In the first and third families the quantum numbers span 
the range n = g, g + 1, . . ., J = g, g + 1, . . . , n and M = 
— J, . . . , J. For the second family n = g + 1, g + 2, . . ., 
J = g + 1, . . . , ri and M = — J, . . . , J. 

The solutions given here are valid in both the north 
pole gauge, Eq. ([3]), and the south pole gauge, Eq. @, if 
we interpret As and Wjm appropriately. 



E. Exceptional modes 

In the charged vector sector there is a new set of modes 
provided g > 1 and these have no counterpart in the neu- 
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tral sector or for g = 1/2. These exceptional modes occur 
for J = g — 1 and M = — J, . . . , + J. The radial quan- 
tum number n becomes continuous as does the frequency 
spectrum. The construction of these modes is described 
in Appendix [B] where we also describe the radial func- 
tions £(%). The explicit form of the modes is: 

a x = 

a g = ^sinxexp[i(M + 5)^(sin6') s+M - 1 (l + cos6l)- M 
a = ^sinxexp[i(M + .g)0](sin6l) 9+M (l + cos6»)- M . 

(30) 

The frequencies of these solutions form a continuum, 
A 2 G (—oo, +oo). 

It may seem surprising that the spectrum is a contin- 
uum on a finite space S 3 . This result is an artifact of 
treating the monopolcs as point objects. As explained 
in Appendix [5] this leads to a singular potential in the 
mode equation for the J = g — 1 angular momentum 
channel; the singular potential in turn is responsible for 
the continuum character of the spectrum. If we worked 
with a model in which the monopolcs had structure the 
singularity would be softened and we would presumably 
obtain a discrete spectrum. However we still expect that 
there would be both bound and unbound modes with A 2 
less than and greater than zero respectively. The precise 
spectrum would depend on the assumed structure of the 
monopole core. Fortunately we do not need to determine 
the exact spectrum for the problems we wish to consider 
in this paper. 



F. Zero modes 

So far we have been looking at modes which are com- 
pletely regular at the poles \ — 0, ir of the S 3 . If we 
admit modes that diverge at the poles, we can construct 
zero modes that can be interpreted as translations of the 
monopoles. 

If we rotate the monopole and antimonopole by a small 
amount while keeping them antipodal the change in the 
field configuration will be a zero mode — a zero frequency 
solution to the linearized Yang-Mills equations. It is a so- 
lution because the original and displaced monopole fields 
are both exact solutions to the full Yang-Mills equations; 
therefore their difference should be a solution to the lin- 
earized equation for small displacement. The frequency 
vanishes because both original and displaced configura- 
tions have the same energy. These zero mode solutions 
are constructed explicitly in Appendix[B] There are three 
independent zero mode solutions corresponding to the 
three orthogonal directions in which we can infinitesi- 
mally displace the monopole- antimonopole pair. Explic- 
itly these solutions are given by: 



First Zero Mode: 



,(3) 

(3) 
>0 

(3) 



o 



cosx 
sinx 



(1 - cos 2 9). 



(31) 



Second Zero Mode 
„(3) - 



cos 9 



x 

'() 

.(3) 



— (- 

sin 2 x ^ sm 
cos x ( cos 9 — cos 2 



sm i 



sm x ^ sin 9 
cos x 1 1 — cos ( 



sm < 



sinx 



— sm o cos t 



cos (f> (32) 



Third Zero Mode: 



,(3) 



,(3) 



,(3) 



— (- 

, 2 vV 



cos 9 



sm x 
cos x / cos ( 



sm ( 



sin 2 9 



sinx 

cos x ( 1 — cos 
sinx 



sin 9 



sin 9 cos 9) sin</> (33) 



The zero modes lie in the uncharged sector because 
they correspond to the difference of two field configura- 
tions with the same charges. Note that the zero modes 
are singular along the line 9 — tt. This is because we have 
worked in the north pole gauge Eq. ([3]) . If we had worked 
in the south pole gauge, Eq. ([J) the solutions would have 
been singular along 9 = 0. The zero modes did not arise 
in the previous analysis of the neutral sector because we 
had previously restricted attention to solutions that were 
globally regular. 

As the three zero modes listed above are associ- 
ated with displacement of the monopole and the anti- 
monopole, they are not present when the background is 
trivial. It is worth noting that there is another solution 
with zero frequency, not associated with displacement, 
and therefore present also for the trivial case: this is the 
neutral scalar sector solution (fl9l) with n = J = M = 0. 



G. Gauge condition 

The gauge conditions, Eqs. (JTSJ) and (jTSJ) , do not mix 
the neutral and the charged sectors but, unlike the equa- 
tions of motion, do mix the scalar and vector compo- 
nents. Using the solutions of the previous section, we now 
find linear combinations of those solutions that satisfy 
the gauge conditions. In other words, we find 4- vectors 
of the form (aa^ , PjS^ (j)), where a and /3j are con- 
stants, and j = 1,2,3 labels the 3 vector solutions we 
have found. The coefficients a and (3j are determined by 
inserting these linear combinations into (|15p or (|16[) . 
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We find three solutions 
(0, aW(2)), and (0 

combination, (iuia 



a 



(fc) 



( ■ Ck) 
However. 



a«(3)) ^ 
( fc ) a( fe )(l)) i s p Ure gauge 



a«(l)), 
the first 



This 



ran 



be easily seen as it is the gradient of a scalar function. 
As expected, there is no longitudinal mode for a massless 
gauge field. Discarding this solution, we are finally left 
with 2 physical (transverse) modes for each sector. 
A similar analysis of the zero modes shows that the 



three solutions, Eqs. (01 



combined with a 



(3) 



0. 



satisfy the gauge condition and are physical. The last 
zero mode, ([TO)) with n = J = M = 0, is pure gauge for 
the monopole-pair background, as well as for the trivial 
background. 



V. SUMMARY OF PHYSICAL MODES 

After eliminating the combinations of eigcnmodes 
that do not satisfy the gauge condition and also the 
combinations that are pure gauge, we are left with the 
following physical eigcnmodes and energies. 

Trivial Background: In the absence of the monopolc- 
antimonopole pair all three directions in group space (± 
and 3) look alike, and each of them has 2 physical eigen- 
modes given by Eqs. ([23)) and ([2"6")) with = (where 
k = ±, 3). Both eigenmodes have the same frequency, 

A = n + 1 , (34) 

where n = 1,2,3,..., J = 1, 2, . . . , n and M = — J, . . . , J. 
There aren't any physical zero modes. 

Monopole-antimonopole background: The neutral 
sector (which is direction (3) in group space) has 3 
zero modes with A = 0, and two physical eigenmodes, 
Eqs. (J23J) and ([2"g]). both with energy 



A 



1 



(35) 



where n = 1, 2, 3, . . ., J = 1, 2, . . . , n and M = — J, . . . , J. 
The (+) and (— ) directions of group space arc similar, 
and each has 2 physical eigenmodes given by Eqs. I|28p 
and ([2"9")> . again with = 0. The energy eigenvalue for 
both modes is 



A = 



1-J + a 



(36) 



where n = g + 1, . . ., J = g + 1, . . . , n and M = 
— J, ...,+J for the physical solution corresponding to 
Eq. ([2"5| . and n = g, g + 1, . . ., J = g, g + 1, . . . , n and 
M = — J, . . . , + J for the physical solution corresponding 
to Eq. ([29)) . The parameter a is defined in Eq. ([21")) . 
For g > 1 there is an additional exceptional mode, given 
by Eq. ([50)) . The eigenfrequencies of these solutions form 
a continuum, A 2 6 (— oo,+oo). 



VI. CLASSICAL STABILITY 

A key consideration is whether the modes are stable. 
The diagnostic for stability is whether the frequency of 
a mode is real or imaginary, or equivalcntly, whether the 
frequency squared is positive or negative. 

We find that the lowest frequency for g = 1/2, which 
occurs in the charged vector sector for n = 1/2 and J = 
1/2, is real and positive: 



min(A) 



V3 + 1 



1.37 



(37) 



(see Eq. ([2"5)) or ([2"9")0 . Hence a monopole with g = 1/2 
is classically stable. 

In the case of g > 1 monopoles there is a continuum of 
exceptional modes, Eq. ([311)) . that have A 2 < and are 
hence unstable. This shows that monopoles with g > 1 
are classically unstable, a result first obtained by Brandt 
and Neri [llf . 

To visualize the instability of g > 1 monopoles, con- 
sider the g = 1 case. Then we have J = g — 1 = and 
hence M = 0, and we find 



fe<f> = dga^ — 8$ 



^sinxe l0 (l-cos6l) (38) 



Therefore the unstable mode develops a radial magnetic 
field component at 9 — it but none at 9 = 0, suggest- 
ing that the instability of g — 1 monopoles is toward a 
spherically asymmetric configuration. 



VII. QUANTUM EFFECTS 

We now turn to the quantum correction to the energy 
of a monopole antimonopole pair. We consider only the 
case g — 1/2 since for higher g values the configura- 
tion is classically unstable. Our objective is to analyze 
the divergent zero-point energy in the monopole back- 
ground. From this divergent contribution we make in- 
ferences about how the energy of the configuration scales 
with the distance R between the monopole-antimonopole 
pair. A faster than linear growth would suggest confine- 
ment. On the other hand if the energy scales more slowly 
or if it decreases, but slower than 1/R, it is still signifi- 
cant because it suggests that at large distance the energy 
of a pair is dominated by quantum effects. 

We aim to compute Eq , the difference in the zero-point 
energy of the monopole-antimonopole and the zero-point 
energy in the absence of the monopole-antimonopole, 



e q r = E A ■ E A 



modes 



modes 



(39) 



trivial 



where the mm subscript refers to A in the background 
of the monopole-antimonopole pair, and the trivial sub- 
script refers to A when there are no monopoles i.e. the 
trivial vacuum. Once we have Eq, we wish to study its 
dependence on the number of summed modes r c , which 
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we also relate to the distance between the monopoles by 
r c = PR where P is an ultraviolet momentum cut-off. 

Ideally wc would sum over modes in the mm sector that 
map onto the modes in the trivial sector. Unfortunately 
it is not clear there is such a mapping. Wc might consider 
bringing the monopole and antimonopole together until 
they annihilate and watch the modes evolve through this 
process. Unfortunately the modes are highly degenerate 
when the monopole and antimonopole are antipodal or 
coincident and so it is impossible to establish a conti- 
nuity between modes at these two extreme points. Al- 
ternatively one could imagine keeping the monopole and 
antimonopole fixed and study how the modes evolve as g 
is turned off. The problem is that for consistency g must 
be half-integer quantized so we cannot continuously vary 
g. Although physically variation of g is impossible, we 
cannot rule out that there is some way to establish a 
mapping by analytically continuing g from 1/2 to zero. 

For the present work, we have subtracted the mode 
energies by arranging the modes (indexed with k) in as- 
cending order in energy, A(fc + 1) > A(fc), and fixing the 
number of summed modes r c , 



E Q R = £A(fc) -J>(A) 



fe=i 



k=l 



(40) 



trivial 



The ordering of the modes according to the energy is 
as follows. In the trivial background the mode energy 
only depends on the n quantum number and we simply 
order the modes according to n. In the mm background, 
the energy of the mode is n + 1 — J + a. It is easy to see 
that the combination — J + a lies in the interval [—1/2, 0) 
and is monotonically increasing with increasing J. This 
implies that the J dependence of the mode energies sim- 
ply gives a small splitting among different states having 
the same n quantum number. Therefore it is sufficient to 
order the modes primarily according to increasing n and 
then, for the same values of n, according to increasing J. 

It is obvious that the charged sector modes are differ- 
ent depending on whether they are computed about the 
trivial background or the background of the monopole 
antimonopole pair. Thus these modes will contribute to 
the Eq which is the difference in the zero point energy 
with the monopole background and the zero point en- 
ergy with a trivial background. A more subtle point is 
that the neutral sector modes are also different about 
the monopole background. In this case there are three 
zero modes that are absent in the fluctuation spectrum 
about the trivial background, so the first mode of the 
neutral sector in the different backgrounds contributes 
differently, 



A( 3 )(l) 



1 



(41) 



= 0, A< 3 )(1) 

trivial 

Thus the neutral modes also contribute to the difference 
in zero point energies Eq that we wish to calculate. 

It is possible to obtain analytical bounds on Eq and 
these are derived in Appendix [D] Here we show our nu- 
merical results for Eq, obtained after summing Eq. (|40[) . 



Neutral sector 


Charged sector 


n 


multiplet 


n 


multiplet 


1 


6 


1/2 


4 


2 


16 


3/2 


12 


3 


30 


5/2 


24 


k 


2k 2 + 4k 


fc + 1/2 


2k 2 + 6fc + 4 



TABLE I: The first few n quantum numbers of the neutral 
(left) and charged (right) sectors, and the number of modes 
associated with that quantum number. 



Fig. [3] shows a plot of EqR versus r c , plotted in solid 
line. The growing oscillatory behavior in Fig. [3] makes 
it hard to interpret the result for Eq. The oscillations 
originate due to the difference in the degeneracy struc- 
ture of the mm and of the trivial case. The total angular 
momentum of the neutral sector is integer, J = 1,2, ..n. 
For the charged sector, the total angular momentum is 
half integer, J = 1/2, 3/2, ..n. Ignoring the fine splitting 
of J values in the charged sector, the nth multiplet has 
J2j 2(2 J + 1) modes. While the form is similar for both 
the neutral and charged sector, the relevant J span is dif- 
ferent, resulting in multiplcts of different sizes. The first 
few n's and their associated multiplcts of both sectors are 
listed in Table HI The result is that each sector reaches 
saturation of an energy level at a different value of T c . 
The energy sum is dominated by the charged sector for 
the monopole-antimonopole background and by the neu- 
tral sector for the trivial background, creating the zig-zag 
effect of Fig.[3J The enveloping cone, plotted with dotted 
lines, is calculated at r c values where either the neutral 
or the charged sectors are saturated and jumping an en- 
ergy level. Fig. [3J also shows the average of Eq (dashed 
line). As each oscillation is almost linear, we average 
only the extremal points, which arc on the enveloping 
cone. Naming the extremal points (Ti,EQi), the average 
is calculated according to 



T(avcragc) 



E, 



(average) 



r. 1 _ 1 + 2r. i + r. 1+1 



(42) 



Eqi-i + 2EQi + Eqi + i J (43) 



In Fig. U we plot average of EqR on a log- log plot and 

1/3 

find (EqR) cx — r c . We have also plotted the envelop- 
ing curves themselves on a log-log plot in Fig. The 
result is that the amplitude of the oscillations grow in 
proportion to ; the lower enveloping curve oc — T 2 / 3 , 
and the upper bound cx +r 2 / 3 . 

A more complete renormalization analysis would re- 
quire us to evaluate other divergent contributions to the 
energy of the configuration and to extract from them a 
finite correction to the energy of the configuration. We 
leave this problem open for later study. Such an anal- 
ysis is needed for a complete evaluation of the leading 
quantum correction to the pair configuration energy. 
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FIG. 3: The energy difference Eq between the mm vacuum 
and the trivial vacuum is found to oscillate as a function of 
the number of modes (r c ) included in the energy sum for 
each sector. Also drawn in are the enveloping curves for the 
oscillations and the mean. 




FIG. 4: Smoothing the oscillatory behavior shown in Fig. [3] 

1/3 

reveals that on average EqR varies as — r c as shown in this 
plot. The smoothing was done as explained in the text. 
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FIG. 5: The amplitude of the oscillations shown in Fig. [3] is 
shown by this logarithmic plot to vary as . 



VIII. CONCLUSIONS AND DISCUSSION 

We have evaluated the complete set of vector spher- 
ical harmonics in the background of a monopolc- 
antimonopole pair on the antipodes of an S 3 . We have 
also solved the radial problem and thus have a complete 
solution to the fluctuation problem. From our results, 
we can confirm the Brandt and Neri instability for g > 1 
non-Abclian monopoles (l4| . The unstable eigenmodes 
are determined explicitly and are found to be asymmet- 
rical in the 9 direction, suggesting an instability which 
might have an interpretation in terms of flux confine- 
ment. 

We show that non-Abelian monopoles with minimum 
charge (g = 1/2) do not have a classical instability. We 
have constructed the complete set of fluctuation modes 
and their eigenvalues, and we use these results to find the 
energy of the zero point fluctuations over and above the 
energy of such fluctuations in the trivial vacuum. The re- 
sulting energy contribution Eq is found to oscillate with 

2/3 

an amplitude that grows as T c where r c is the number 
of modes that we include in the sum over fluctuations. 

1 /3 

The mean value of Eq however is proportional to — IV . 
If the number of modes in the sum is limited by an ultra- 
violet momentum cutoff, P, then r c = PR and we find 
that E Q oc -i?~ 2 / 3 . In other words, the mean contribu- 
tion of the vacuum fluctuations increases with R but goes 
to zero in the R — > oo limit. For confinement we would 
expect Eq ~ R in the R — ► oo limit and so we conclude 
that our results do not provide evidence for quantum con- 
finement of non-Abelian magnetic monopoles, at least on 
S 3 . 

Our analysis does indicate that quantum effects be- 
come important at large R because the zero point energy 
oc — i?~ 2 / 3 goes to zero more slowly than the Coulomb 
energy oc — 1/R. At some critical value of R, quantum 
effects overtake classical effects. However, such a con- 
clusion relies on the assumption that a full treatment 
of the renormalization of the monopole mass does not 
cancel out the — i?~ 2 / 3 dependence. In order to check 
this assumption, we would need to quantize regular non- 
Abelian monopoles e.g. 5/7(3) monopoles as discussed 
in the introduction. 

A puzzling feature of our analysis is the occurrence of 
oscillations in Eq as a function of the mode cutoff T c . 
We are convinced that these are present in our scheme 
of comparing energies in the mm and trivial sectors by 
taking an equal number, r c , of lowest lying states in ci- 
ther sector. The underlying reason is that the degener- 
acy structure of the modes is different with and without 
monopoles. Mode degeneracies depend on the symme- 
tries of the system and it is clear that the mm and trivial 
sectors have different symmetries under rotations. So it is 
not surprising that the degeneracies in the two cases are 
different. Once we accept the existence of different de- 
generacies, and we compare sums with the same number 
of modes, then the steps in the principal quantum num- 
ber, n, in the sums over modes, arise at different values 
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of r c for the cases with and without the monopoles. For 
some T c the sums with the monopoles give a larger result 
and for other r c the sum in the trivial vacuum is larger. 
While we understand the oscillations mathematically, we 
cannot rule out the possibility that if we were to track 
the flow of modes in transitioning from one sector to the 
other, say by bringing the monopoles closer and letting 
them annihilate, it may lead to a different subtraction 
scheme. We hope that future work will shed more light 
on these issues. 

Finally we wish to remark that the quantum in- 
teraction of non-Abelian magnetic monopoles may be 
amenable to analysis by lattice methods and this would 
be an alternative approach to solving the problem. 



vector equations of motion. 



V^ a = d x a x + 2^a x 



1 



sin 2 x 



smx 

cos 9 
Ogag + ^—^a e 



1 



sm 9 sin 9 
A (3) = dl + 2^9. 



d<f,a d 



(Al) 



A sm x 



-A- U + C -^d e + -^-dl] (A2) 
sin x\ sm ^ sin 8 I 
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(A< 3 >a) x = AW 



2 \ _ 9 cosx 
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costf 
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V sm x sm " J 
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smx 
2 



cosf 



sin 2 x sin 2 9 sin ( 



(A3) 
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(A^a) = A* 3 ^ + 2^ d^a x - 8 x a 4 
sm x \ 



2 cos( 
sin 2 x sm ' 



(A5) 



APPENDIX A: DIFFERENTIAL OPERATORS 



In this appendix, we give the explicit form of Eqs. (jll[) - 

(USD. 



1. Neutral sector operators 



Below are the equations of the uncharged sector. The 
(3) index of the vector components a,; is suppressed. The 
operator given in Eq. (|A1[) is used in the gauge condition, 
and the operator in equation (|A2[) is used in the equa- 
tion of motion for the scalar component. Eqs. (|A3|) - (|A5|) 
are the components of the vector Laplacian, used in the 



2. Charged sector operators 



In the charged sector, derivatives are gauged, and 
Zeeman-like terms appear. We use the notation of Aa = 
(1 — cos9), g = eq/in, and A = — igA<p. Again, the 
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(+) index of the vector components aj is suppressed. 



V (+) a= d x a x + 2 



cos x 

( 

sinx 



-a. 



1 



sin 2 x 



do 



cos ( 



1 



sin^" 17 ' sin 2 # 

A (+) = Ql + 2 ^d x + 
K sin x 



sin 2 x 



(A (+) a) x = A 



(+) 



sin 6* sin" 6* 



2 \ o cosx 
a v — 2- 



(A6) 



(A7) 



2 i X ^ • 3 

sin x / sm X 



cos \ i 

09 + -t— pr as + . 2 „ D <t> a <t> 
sm.V ] sin 



(A8) 




sin 2 x sin 2 \ sm ( 



£></> - igdeAtf, (A9) 



(Aa(+)) = A(+>a + 2^ A^x - d x a 4 
sinx \ 



sin 2 x 



COS! 

sin t 



— D^ag 



(A10) 



Note that in the case that g = the equations of the 
charged sector coincide with those of the neutral sector. 



This is Gegcnbauer's equation which has solutions 

regular at i = ±1 only for the quantized values 

A 2 = n{n + 2). These solutions arc the Gegen- 
j + 1 

bauer polynomials G n _j(x) where n — J must be a 
non-negative integer. Thus we obtain the solutions 
£n-J,j(x) Y JM(Q,<f>)exp(iut) where £„-j,j is given by 
Eq. (fTS]) in the text. The quantum numbers are restricted 
to n = 0, 1, . . . and J = 0, 1, ... n and M = — J, . . . , J. 

As noted above the boundary conditions we impose 
are that the solution should be regular at x = and 
X = 7T. The physical basis for these boundary conditions 
is that in the vicinity of the monopole, x = 0, the back- 
ground field strength diverges as (F^ F^) 1 / 2 oc sin -2 x, 
and there is a similar divergence in the vicinity of the 
antimonopole, \ — 7r - I n order for the perturbative ex- 
pansion used in Eq. ([5]) to be valid, the divergence of the 
field perturbation /„„ cannot be stronger than the back- 
ground divergence. If the divergence is weaker than the 
background then 



lim 



fill* f 

J J liv 







x^yF^F^ ^ 
or if it is as strong as the background solution, then 



(B4) 



//// 



lim , 



Const 



(B5) 



In practice we find that the regular solutions to Gegen- 
bauer's equation satisfy these boundary conditions; the 
singular ones do not. In fact the regular solutions all 
prove to be less singular than the background and so sat- 
isfy the more stringent boundary condition (|B4[) . The 
boundary conditions we impose in the other sectors are 
similarly motivated. 

For the trivial background which doesn't have F^ v , we 
demand regularity at x = and x = it. 



APPENDIX B: DERIVATION OF EIGENMODES 



2. Charged scalar sector 



1. Neutral scalar sector 

We wish to find the eigenmodes of the scalar wave 
equation on S 3 . The obvious ansatz 



4 3) =Z(x)Yj M (6,<t>)e iwt 
leads to the eigenvalue equation 



dz 



,cosx 
sinx 



V 



A - 



J(J + 1) 



sin 2 x 



(Bl) 



K = 0. (B2) 



Making the substitutions x — cosx and £ = (1 — x 2 ) J / 2 r) 
leads to 



In this sector we wish to solve the scalar wave equation 
with minimal coupling to the background vector poten- 
tial of the monopolc-antimonopole pair. The solution 
closely parallels the neutral case. The natural ansatz is 
to replace spherical harmonics with Wu-Yang monopole 
spherical harmonics: 



(B6) 



This leads to the same equation for £ as in the neutral 
case, Eq. (|B2[) . but with the replacement J( J + 1) — > 
J(J + 1) — g 2 . To bring out the parallel it is convenient 



to define a as in Eq. (|21[) of the text, so that a(a + 1) = 
J(J + 1) — g 2 and a > 0. Next we substitute x = cosx 
and £ = (1 — x 2 ) a l 2 r\ to find that r) obeys Gegenbauer's 
Eq. (|B3|) with the replacement J — > a. Solutions that 
are regular at both x = ±1 can only be found for the 
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quantized values uj = [n — J + a) (n — J + a. + 2). These 

solutions are Gegenbaucr polynomials G n _ j where n — J 
is a non-negative integer. 

In summary the solutions in the charged sector are 
£n-j.a(x)WjM {9, 4>) exp(ia>t) with frequencies uj 2 = (n — 
J+a)(n—J+a+2). The quantum numbers are restricted 
to the range n = g, g + 1, g + 2, . . . and J = g, g+ 1, . . . , n 
and M = —J,..., J. The condition on n arises from 
the Gegenbauer quantization condition; the limits on J 
and M are the usual ones in the theory of monopolc 
harmonics. 



3. Neutral vector sector 

We wish to determine the eigenmodes of the vector 
wave equation on S 3 . It is helpful to first solve the 
corresponding problem on R 3 . It is natural to seek so- 
lutions of the form f{r)X%(6,<l>), /(r)X^(0, 0), and 
/(r)Xj^(#, (j)). The three families correspond to eigen- 
modes of different polarization. Here Xjm(0>^>) are the 
vector spherical harmonics (22| . They are eigenfunctions 
of total angular momentum (orbital plus spin) with quan- 
tum numbers J and M, They are also eigenfunctions 
of orbital angular momentum squared with eigenvalues 
J(J + 1), (J + 1)(J + 2) and ( J - 1) J respectively. The 
quantum numbers span the range J — 0,1,2,... and 
M = — J, . . . , + J. The case J = is special in that there 
is only one vector spherical harmonic Xqq ; the other two 
X^ and Xq ' vanish. 

By use of the wave equation we determine that the ra- 
dial functions f(r) are spherical Bessel functions of order 
J,J+l and J— 1 respectively. Thus in polar co-ordinates 
the solutions designated /(r)Xy2f(0, </>) are given by 

V r = 

ve = rjj(wr)^-rd < i > Yj M 
sin u 

v,/, = -rjj(ur)sm6dgY JM ; (B7) 
the solutions /(r)Xjy by 

v r = -(J + l)jj + i(ur)Yj M 

Ve = rjj +1 (ur)deYj M 

Vet, = rj J+1 (u>r)d4,YjM\ (B8) 

and the solutions designated f(r)Xj^ by 

V r = j.j-i{u)r)JYj M 

ve = rjj-i(ujr)dgY JM 

v<p = rjj_i(u!r)d<pYjM- (B9) 

The frequency uj 2 is continuous and restricted to the 
range uj 2 > 0. 

For the subsequent generalization to S 3 it is necessary 
to introduce a different fundamental set of solutions. In 



the new set, the first family of eigenmodes is obtained 
by taking the gradient of the modes of the scalar Lapla- 
cian on R 3 , namely jj(ur)YjM- For the second fam- 
ily we adopt the transverse solutions /(r)Xj°^ / given in 
Eq. (|B7[) . For the third family we take the curl of the 
second group of transverse solutions. Explicitly, then the 
first family of solutions are 

v r = d r jj{ur)Yj M 
v g = jj(ujr)d g YjM 
V4, = i.j{ur)d,pYjM- (BIO) 

The second family is given by Eq. (|B7j) . The third family 
of solutions, obtained by taking the curl of the second 
family, are 

v r = -jj(ur)J(J + 1)Y JM 
r 

vg = d r [rjj(ujr)}deYjM 

v<t> = d r [rjj(u}r)]d4,YjM- (Bll) 

In the special case J = the transverse solutions and 
their curls vanish; only the gradient solutions survive. 
This is consistent with the expectation that there is only 
one polarization in this exceptional angular momentum 
channel. 

In summary the first independent set of solutions is 
given by E qs. ijBTjl . (|B8j) an d (|B9|) ; the second set by 
Eqs. |[E10]) . p37|) and (|B"TTj) . By making use of re- 
cursion relations for Bessel functions we can show that 
the gradient solutions Eq. (|B10[) are the superposition 

/(r)X^ + /(r)X^; the curl solutions Eq. (|BTT|) are 
the superposition - J/(r)X^ + ( J + l)/(r)X^. Thus 
the two alternative sets of eigenmodes are seen to be 
equivalent. 

We are now ready to tackle the problem on S 3 . For the 
first set of eigenmodes we try the gradient of the scalar 
modes leading to the ansatz 

4 3) = d xUjY.JM 

af = tnjd+Yju (B12) 

By applying the S 3 vector Laplacian to these modes we 
determine that their eigenvalues arc n(n + 2). 

For the second set of solutions we seek purely trans- 
verse solutions. By analogy to Eq. (|B7|) we make the 
ansatz 

4 3) = o 

af = -f(x)em08oYj M (B13) 

Application of the S 3 vector Laplacian shows that these 
functions are eigenmodes if we take f(x) ~ * s i n X?nj(x)- 
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These modes are found to have frequency oj 2 = n(n + 
2) + l. 

The third set of solutions are obtained by taking the 
curl of the second transverse set. Again application of 
the S 3 vector Laplacian shows that the resulting modes 
are eigenfunctions with frequency lu 2 = n(n + 2) + 1. 

By analogy to R 3 we see that for the second and third 
set of solutions n = 1,2,3,... and J = 1,2,3,... and 
M = —J,..., J. For the first set, n — 1,2,3,... and 
J = 0, 1, 2, ... and M = -J, ... , J. 

This completes the derivation of the modes enumer- 
ated in section IV. Note it is possible to systematically 
rederive these eigenmodes for the vector Laplacian on S 3 
using the theory of so(4) representations just as in the 
scalar case. Although we will derive the relevant symme- 
try generators below we will not carry out this analysis 
here. 



4. Charged vector sector 



As noted in Sec. IIII| the eigenmode problem in the 
charged vector sector may be interpreted as the non- 
relativistic Schrodinger equation for a spin one particle 
confined to a 3-sphere and moving in the magnetic field of 
a monopolc-antimonopole pair placed at opposite poles 
of the 3-sphere. The particle is minimally coupled to the 
magnetic field due to it's charge as well as Zeeman cou- 
pled to the magnetic field due to its magnetic moment. 

Again, as a prelude, let us solve the simpler problem 
of such a particle moving in R 3 in the field of a sin- 
gle monopole at the orig in, a problem first investigated 
by Brandt and Neri [ljj. In ref [l4| the focus was on 
finding unstable modes by looking for imaginary eigen- 
frcquencics. We extend that work by enumerating all 
stable modes and deriving explicit expressions for the 
eigenfunctions. Our primary interest in the R 3 problem 
is that it provides an important testing ground for the 
ansatze that we will later deploy on S 3 . 

In the absence of the Zeeman term it is natural to seek 
three families of solutions of the form /(r)X 



(o) 



/HXS(W), and /(r)X&(M). Here Xj m{0,4>) are 
the monopole vector spherical harmonics [lg|. They are 
eigenfunctions of the total angular momentum (orbital 
plus spin) with quantum numbers J and M. They are 
also eigenfunctions of the total orbital angular momen- 
tum with eigenvalue J(J+1), (J+l)(J + 2) and (J-l)J 
respectively. For g > 1 the quantum numbers span the 
range J = g — 1, g, g + 1, . . . and M = — J, . . . , J. How- 
ever for J = g — 1 there is only one vector spherical 
harmonic X^ + ); for J = g there are two polarizations 
X<+) and X(°). For all higher values, J = g + 1, g + 2, . . ., 
all three polarizations exist. Similarly for g = 1/2 the 
quantum numbers span the range J = 1/2, 3/2, . . . and 
AI = —J,..., J. However for J = 1/2 there are only 
two polarizations X' + ) and X^ ). For all higher values, 
J = 3/2, 5/2, . . ., all three polarizations exist. 

The Zeeman term couples spin and orbital motion but 



is invariant under total angular momentum. By the 
Wigner-Eckart theorem therefore it can only couple total 
angular momentum multiplets with the same quantum 
numbers. Thus we conclude that for a given J and M 
there will still be three families of solutions but they will 
be superpositions of the form 



f(T) 



a X 



(0) 
JM 



b X 



(+) 

JM 



c X 



(-) 

JM 



(B14) 



The coefficients a,b,c are constants not only in that they 
are independent of (r, 0, <fr) but also of M (in accordance 
with the Wigner-Eckart theorem). They are computed 
by solving for the eigenvectors of a 3 x 3 matrix that is 
determined by substituting the ansatz, Eq. (|B14[) . into 
the eigenvalue equation for the charged vector sector. 
This matrix is given explicitly by Brandt and Neri; it 
can be efficiently derived making use of identities given in 
Ref. Substitution of the ansatz, Eq. (|B14[) . into the 
eigenvalue equation also shows that f(r) satisfies the ra- 
dial equation for a spinless free particle in non-relativistic 
quantum mechanics, 



0. 



2 1(1 + 1) 



A 2 )/ 



0. 



(B15) 



with the modification that the angular momentum I — > 
a — 1 or a or a + 1 corresponding to the three families 
of solutions. For definiteness we shall call these families 
s, t and u respectively. Here a is given by Eq. (|2ip of the 
text. 

Explicitly the solutions are given by 



V r = [JiJ+^-g^j^W 

v g = (a + l)rJ a _ 1 d e W 
= (a + l)rj a _ 1 D c i > W 



(B16) 



for the s-polarization; 

v r = 



ve = r Ja \[J(J+l)-g 2 ]—D^W + gdeW 
v^ = rj a (gD^W - [J [J + 1) - g 2 \i sin 6d e W) 



(B17) 



for the f-polarization; and 

V r = -[J(J+l)-9 2 ]ja+lW 

ve = rj a+1 d e W 
V4, = rj a+1 D<pW 



(B18) 



for the u-polarization. Here j a (Ar) arc spherical Bcssel 
functions of order a. For brevity we have written the Wu- 
Yang harm onics W JM {0,<l>) as W in Eqs. (|B16j) . (|B17|) 
and (|B18p . These modes form a continuum with fre- 
quency A 2 > 0. These three families of solutions are 
analogous respectively to to Eqs. (|B9|) . (|B7|) and (|B8|) to 
which they reduce for g = 0. The t- modes derive their 
name from being transverse. 



14 



There are two complications we have overlooked in 
Eqs. ([ET7]) . (jBT6l) and (|BT8)) . First we have ignored the 
special values of J for which there arc fewer than three 
polarizations. These special cases can be analyzed sim- 
ilarly. For example for case that J = g we must drop 
the c term in the ansatz since there is no vector spherical 
harmonic X^ - ) in that angular momentum channel. The 
upshot is that there are only two families of solutions 
in this angular momentum channel. These turn out to 
be the s and u solutions given above. It can be shown 
that the t solution is identically zero for J = g. For the 
case J = g — 1 we must drop the b and c in the ansatz 
Eq. (|B14[) and there is only one family of solutions. These 
solutions have the form 

v r = 

v g = rfe^q>{i[M + g}(f>){sin9) 9+M - 1 (l + cos6)- M 

w = ir f cxp{i[M + g ](t)) (sin ey j+M (i + cos ey M . 

(B19) 

Here M = —(<?— 1), ...,(</— 1) and f(r) obeys E^JHTB 
with 1(1 + 1) = — g. The angular factors in Eq. (|B19|) are 
taken from the angular dependence of the vector spheri- 
cal harmonics Xg^\ AI (9, (/)). 

The second complication is that in non-relativistic 
quantum mechanics 1(1 + 1) is a non-negative integer and 
the term involving I in Eq. (|B15[) may be interpreted as a 
centrifugal barrier. Mathematically this means there are 
no bound state solutions; the continuum states are spher- 
ical Bcssel functions. Here since we replace I — > a— 1 or a 
or a + 1 it is a concern that 1(1 + 1) may become negative 
depending on the particular value of a. For the negative 
case Eq. (|B15[) would have the interpretation of a parti- 
cle in an attractive inverse cube central force potential. 
According to standard lore [23| in this case for a weakly 
attractive potential, one for which > 1(1 + 1) > —1/4, 
there are no bound states and the unbound continuum 
states are still simply spherical Bessel functions. For 
1(1+1) < —1/4 there is a continuum of bound states 
with A 2 < as well as a continuum of unbound states 
with A 2 > 0. 

The continuum of bound states is an artifact of treat- 
ing the monopole as a point. If we worked with a model 
in which the monopole had structure the singularity in 
the attractive inverse square potential would be softened 
at the origin. Presumably this would lead to a discrete 
bound state spectrum which is generally expected for 
non-singular potentials. 

For J = g, we see that 1(1 + 1) can indeed be negative 
for the s polarized states for g = 1/2, 1 and 3/2. However 
in all these cases it is not so negative as to form bound 
states (recall that for an s polarized state I — > a — 1 
where a is given by Eq. (|21j) . It is now easy to verify 
that -1/4 < 1(1 + 1) < for < g < 2; otherwise 1(1 + 
1) is positive). This leads to the important conclusion 
that g — 1/2 monopoles are stable (recall bound states 
correspond to instability). For J = g — 1 however 1(1 + 
1) — > — g allowing the formation of bound states. Thus 



monopoles with g > 1 exhibit instability in the J = g — 1 
channel. 

In summary Eqs. (|B16[) . (|B17|) and (|B18|) are the eigen- 
modes needed to analyze the stability of a monopole in 
R 3 . For g = 1/2 they are stable and represent a complete 
enumeration of modes. For g > 1 there is an additional 
branch of unstable modes with J = g — 1 . 

The chief virtue of the derivation above is that it is 
systematic and can be counted upon to yield a com- 
plete set of solutions. As in the neutral vector sector 
we now rederive the eigenmodes in a second more in- 
tuitive manner. The merit of this second derivation is 
that it produces slightly simpler expressions and gener- 
alizes readily to S 3 . This second set of eigenmodes is 
also organized into three families. The first family is 
obtained by starting with the gradient solutions of the 
neutral vector sector, Eq. (|B10j) and making the mini- 
mal coupling substitution 3$ — > D^, replacing the spher- 
ical harmonics with Wu-Yang harmonics and substituting 
J(J + 1) — ► J(J + 1) — g 2 . Finally we need to replace the 
Bessel functions jj(ujr) with j a (ujr). The third family is 
obtained by making the same replacements to the curl so- 
lutions of the neutral vector sector, Eq. (|B11|) . Curiously 
the same substitutions made to the transverse solution of 
the neutral sector do not lead to a solution in the charged 
sector. For the second family we therefore retain the t 
modes of Eq. (|BT7|) . 

We are now ready to tackle the problem on S 3 . For the 
first family of eigenmodes we start with the gradient so- 
lution of the neutral vector sector Eq. (|B12|) and replace 
spherical harmonics with Wu-Yang harmonics and make 
the minimal coupling substitution — > D^. This is 
found to be a solution in the charged sector provided we 
also modify the radial function £ n -J,j(x) ~> £n-J,a(x)- 
The eigenvalues are found to be (to— J+a)(n— J+a+2). 
For the third family we work with the curl solution of 
the neutral vector sector Eq. (|26|) and in addition to the 
same substitutions we replace J(J + 1) — > J(J + 1) — g 2 . 
This procedure too yields a solution with eigenvalue 
A 2 = (to + 1 — J + a) 2 . Obtaining the second family 
of solutions is more challenging. We give an argument 
below that there have to be transverse eigenmodes. Bol- 
stered by this argument we make the ansatz 

4+> = o 

4 +) = [(J(J+l)- g 2 ]^(D^W)f + gd 6 Wf 
sm# 

a { + ] = gD^Wf - i[J(J + 1) - g 2 ] sm9d$Wf (B20) 

/ is a function of x, W is the Wu-Yang harmonic 
Wjm(&, </>)■ This ansatz is motivated by the form of 
transverse modes in R 3 and is found to yield solu- 
tions provided / = sinx£,n-j,a(x)- The eigenvalues are 
(n + 1 - J+ a) 2 . 

As in R 3 it is possible to show that the second trans- 
verse family does not exist for J = g. There is also 
a single branch of unstable solutions for J = g — 1 for 
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g > 1. These solutions have the form 
a x = 

a g = sm X £exv[i(M + g)(/)}{sin9) 9+M - 1 (l 
a4 , = isinx£exp[i(M + g)<£](sin0) s+M (l 4 



cos 0) 



COS ( 



-Af 



where the radial function £(%) satisfies Eq. (|B2|) with 
J(J + 1) — > —p. A simple analysis reveals that both so- 
lutions to this equation diverge as £ oc l/y^X as X ~~ * 
and similarly as x ^ t. Since both solutions behave ac- 
ceptably at the boundary points there is no quantization 
or restriction on the frequencies A 2 . Hence we conclude 
that in the J = g — 1 channel there is a continuum of 
bound states with A 2 < and a continuum of unbound 
states with A 2 > 0. 

It may seem counterintuitive that the spectrum of 
modes is continuous on a finite space. Similar to the R 3 
case, this result is an artifact of treating the monopoles 
as points. If we had worked in a model in which the 
monopole cores had structure the singularity of the po- 
tential in the radial equation would be softened at the 
north and south pole of S 3 . This would lead to a discrete 
spectrum, which is generally expected for non-singular 
potentials. However we would still expect both bound 
and unbound modes. 

This concludes the derivation of the charged sector so- 
lutions and the range of allowed quantum numbers that 
are given in section IV. 

The main difficulty with solving Eqs. (|A8|) . (|A9|1 and 
(|AT0|) is that they are a system of coupled partial differ- 
ential equations. A puzzling aspect of our solution above 
is that we were able to obtain two solutions by making 
simple substitutions in the neutral vector solutions but 
the transverse solution did not yield to this strategy. For 
this reason we would like to present another line of ar- 
gument that demonstrates there must be a transverse 
solution. This approach can also be applied to the neu- 
tral sector but here we concentrate on the more vexing 
charged case. 

We begin with the gauge condition (|16[) . re- written as 



dgag 



cos 9 



+ 



1 



sin 2 x Rodta t - d x a x - 2 



cosx 

( 

smx 



(B21) 



where we have suppressed the (4-) index. By substi- 
tuting this in the x component of the vector Laplacian, 
Eq. (|A8[) . we eliminate the dependence on ag and ou, 



-4 a, 



sm% 



sm x 



\- ( d e + + -K-Dl ) «x = ^ t a t .(B22) 

sm f sm 



sm x 



Following the physical analogy to the Schroedingcr equa- 
tion, we now look for 3 solutions to the above equation. 



Repeatedly, our aim will be to somehow force separation 
of the x an d the (9, (j>) dependence. We will assume a 
trivial time dependence of the form exp[iwt]. Given that 
time dependence and for a general ui, the above equa- 
tion allows one of three options - a) both at and a x are 
non-zero, b) at is zero but a x isn't, or c) at and a x are 
zero. 

A reasonable first ansatz is at = Ca Xl C being a 
constant. We posit a separation of variables a x = 
P(9, 4>)f{x) m Eq. (|B22[) . solve this equation for a x , then 
use the resulting solution in Eqs. (|A9[) and (|A10[) to de- 
termine ag and a$ yielding the first solution Eq. (|27[) . 

Next we attempt a solution with at = 0. Eq. (|B22[) 
then gives a solution for a x . Armed with this solution 
we return to Eqs. (|A9[) and (j A10|) and posit that ag and 
a^ have the same x dependence. To be specific we try 



(10 



kP(0,(j> 
G{9,4>) 



V fiX) 

sm x 
dixV" 



sin 9 

a = isin9H{9,4>)g{x)e iuJt (B23) 

where k is a constant. Plugging Eq. (|B23[) into the gauge 
condition gives the relations that g = d x f and k sm 6P + 
dgG 4 iDfkH = 0. Using these relations in Eqs. (|A9|) and 
(|A10j) gives the third solution Eq. (|2"5|) . 
Finally let us attempt a transverse ansatz 



= 
_ G(9, 



smt 

a = ism9H(9,<j))g(x)e tuJt . (B24) 

Plugging this into the gauge condition, Eq. (|B21[) gives 
the relation that dgG = iD^H . We assume that D^H = 
i(M +gA)H. Substitution of these results into Eqs. (|A9|) 
and ([ATO]) yields 



d 2 x + A 2 
cos 9 



fix) 



25- 



sm x 
doA 



(B25) 



sin 9 M+gA 
[M + gAf 



+N- 



G = 



(B26) 



where we used N = J (J + 1) — g 2 for brevity. While 
Eq. (|B25|) is easily transformed into the Gegenbauer 
equation, the angular Eq. (|B26j) is less readily solvable. 
Nonetheless because the variables separate, these two 
equations demonstrate the consistency of the transverse 
ansatz Eq. (|B24[) . Eq. (|B26|) is essentially the Wu-Yang 
monopole harmonic equation but with an extra term 
2g(dgA)/(M + gA)dg that makes the equation difficult 
to solve. However we can a posteriori verify that the 
solution Eq. (|28|) is not only consistent with the trans- 
verse ansatz Eq. (|B24[) but also an explicit solution to 
Eqs. (|B25|> and jB26}. 
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5. Zero-Modes 

Here we use the group theoretic technique discussed in 
Appcndix[C]to find the translational zero mode solutions. 

First let us develop expressions for rotation generators 
acting on vector fields in R 4 . If we initially represent vec- 
tors in terms of their cartesian components («i, V2, U3, W4) 
then the generator of rotations in the X3-X4 plane is given 
by 



i cos 9d v — i sin 9 



cos x 
s'mx 



do 



000 
000 
000 

i 






—i 




(B27) 



This result follows from Eqs. (|C1[) and (|C6|) . We would 
rather work with polar components which are related to 
cartesian components via 



(B28) 







v r 


V2 


= u 


v x 


V3 




ve 






. v 4> . 



where the transformation matrix U is given by 



U 



-scj)/rsxs9 
C(j)/rsxsO 



sxs9ccj> cx&9c(j)/r c9c<p/rsx 
sxs9s(j> cxs9s<fi/r c9s(j>/rsx 
sxc9 cxcQ/r —s9/rsx 
cx sx/r 

(B29) 

Here for brevity we have written sin 9 = s9, cos x = C X 
etc. In terms of polar components N3 is given by 
U^ 1 N3 (cartesian) U . A simple calculation reveals the po- 
lar form 



N* 



i cos #<9 V — i sin 9 



cosx. 






— isinf 




smx 



1 sin 9 1 sin 2 x 
— i cos 9 cosx 




(B30) 



We are interested in vectors that live on the unit sphere in 
R 4 . Such vectors have no radial component and we may 
therefore discard the first row and column of the matrix 
in Eq. (|B30[) . The 3x3 matrix differential operators 
corresponding to rotations in the other planes may be 
constructed similarly. 

Now let us construct the zero modes. The background 



monopole field has polar components A x 



A s = 



and A,p = g(l — cos 9). The zero modes are obtained 
by considering the change in this field under infinitesi- 
mal rotation in the X1-X4, X2-X4 and X3-X4 planes, or, 
in other words, by application of Nx,N% and N3 to the 



background field configuration. Application of N3, for 
example, yields (N 3 A) X = (N 3 A) g = and (A 3 A)^ = 
—ig(cosx/ sin%) sin 2 9 which is the first zero mode listed 
in section IV. The other two are obtained by application 
of Ni and N 2 . 



APPENDIX C: GROUP THEORETIC 
DERIVATION OF MODES IN NEUTRAL 
SCALAR SECTOR 



It is instructive to rederive this result using the rota- 
tional symmetry of the problem. Readers who are not 
interested in a group theoretic rederivation may skip this 
section but should still examine Eqs. (|C1[) and (|C6[) and 
read the three sentences preceding Eq. (|C1[) and the en- 
tire paragraph containing Eq. (|C6[) as these results will 
be used to construct the zero mode solutions. 

Picture the sphere S 3 embedded in a four-dimensional 
space. It is easy to write down the 4x4 matrices corre- 
sponding to rotation generators in the x\ — X2, X2 — x 3 , 
X3 — xi and x\ — x&, X2 — X4, x 3 — x± planes. We denote 
these generators Jx, J2, J3 and Ni,N2,N 3 respectively. 
For example 



■h 



-i 

1 





N3 = 







-! 

i 



(CI) 



J3 and N3 commute and together constitute a Cartan 
subalgebra. By computing the commutators of all six 
matrices we can obtain the structure constants of the 
so(4) algebra. It is convenient to define 



K 4 



(•A 



N+); M+ 



A-_ = -(J_ + JV_); M_ = - 
K3 = ^(J 3 +N 3 ); M 3 = i(J 3 



(J+ - N+) 
(J--N-) 
N 3 ) 



(C2) 



where, as usual, ,7+ = J\ + i J 2 and J_ = J\ — i J 2 and 
the same for N. By using the explicit matrices it is easy 
to verify that K and M separately obey the angular mo- 
mentum algebra and commute with each other. In other 
words so (4) = su(2) © su(2). 

Working with the K,M generators, Eq. (|C2[) . we see 
from the theory of angular momentum that representa- 
tions of the so(4) algebra can be labeled by l\ and I2 
which may take any value from 0, 5, 1, |, . . . The repre- 
sentation (li, I2) has dimensionality (2Zi + l)(2Z2 + l). The 
basis vectors in this representation are denoted \mi,vri2) 
where mi = — Zi, ...,Zi and rri2 = — Z2, ...Z2. The ef- 
fect of the basic algebra elements in this representation 
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is given by 
K 3 \m 1 ,m 2 ) = m 1 \m 1 ,m 2 ) 

K + \m ll m 2 ) = [(Zi + mi + l)(ii — mi)]5|mi + 1,TO2) 

K-\mi,m 2 ) = [(Zi +mi)(Zi — mi + l)]'|mi — 
M 3 |mi,m 2 ) = m 2 |mi,m 2 ) 

M+|mi,m 2 ) = [(l 2 +m 2 + l)(l 2 -m 2 )]i\mi,m 2 + l) 
M-\mi,m 2 ) = [(h + m 2 )(h - m 2 + l)]^|mi, m 2 - 1) 

(C3) 

Instead of the |mi,m2) states it is sometimes conve- 
nient to work with the \JM) basis defined by 



\JM) 



E 



E 



Ct£ m ^hM)\™i,™v) 



(C4) 



mi =— 1\ m2=- 1<2 



where J = \l\ — l 2 \, . ■ ■ , h + l 2 and M = — J, . . . , J. 
Cmfm 2 (h>h) are the Clebsch-Gordan coefficients. There 
are (2Zi + 1) (2Z 2 + 1) of these states, as there should be, 
since the | JM) states are simply an alternative basis for 
the {h,h) representation of the so(4) algebra. In this 
basis the matrices for J = K + M are simple but not for 
N or K or M. 

The natural Casimir invariant for the so(4) algebra is 



C 



N 2 



2K 2 



2M 2 



(C5) 



All the states in an (Zi, Z 2 ) representation are eigenstates 
of C with eigenvalue 2Zi(Zi + 1) + 2Z 2 (Z 2 + 1). 

Square representations where l\ = l 2 = I are of par- 
ticular interest. It is helpful to define n — 21. Thus 
a square representations is labeled by a single integer, 
n, its order. A square representation has dimensionality 
(2Z + l) 2 = (n + l) 2 . For square representations we see 
that in the \JM) basis the quantum numbers span the 
range J — 0, ... ,n and M = — J, . . . , + J. Moreover the 
Casimir invariant is given by n(n + 2). 

Now let us consider wavefunctions ip(9,(f>,x) on S 3 . 
The Hilbert space of these wavefunctions constitutes a 
reducible infinite dimensional representation of the so(4) 
algebra. By considering the rotation of a wavefunction 
we can deduce the generators of rotations in the six fun- 
damental planes. For example 



J3 = — ids, = icos6*<9 v 



i smfl Og 

sinx 



(C6) 



These differential operators obey the so (4) algebra de- 
fined by their 4x4 matrix counterparts. In deriving the 
differential operators it is helpful to recall that the rela- 
tion between the polar co-ordinates of a point on S 3 and 
its four-dimensional cartesian co-ordinates is 



•'■'2 
•i'3 
•i'4 



sin x sm " cos < 
sin x sin 6 sin q 
sin x cos 9 
cos x 



and that J3 = —ix\d 2 + ix 2 d\, A3 = —ix^di + ix^d?,. 

A simple calculation reveals that the Casimir differen- 
tial operator 



C 



sm x 



<V,7 



1 



• 2 • 2 

sm x sm 



-A 



2^5, 

smx 



cos ( 



1 



sin 9 



sm x 



(C8) 



coincides with the Laplacian on S 3 . Furthermore it is 
possible to verify that the eigenfunctions of the Lapla- 
cian £nj(x)YjM(9,4>) that we obtained by separation of 
variables have the following group theoretic interpreta- 
tion: For a fixed n, they constitute the \ JM) basis for a 
square representation of order n. 

That only square representations are realized can be 
demonstrated by seeking functions that satisfy K+ip = 
and M + ip = 0. These are first-order equations and the 
solutions are readily found to be 



ip n = (sin x)™ (sin 0)™ exp(m</>) 



(C9) 



where n is a non-negative integer. Since these are the 
highest weight states of the representation we can de- 
duce the (Zi,Z 2 ) values of the representation by applica- 
tion of K3 and M3 to ipn- By explicit calculation we 
find K^ipn = M^ipn = (n/2)ijj n revealing that the rep- 
resentations under consideration are indeed square with 
Zi = Z 2 = n/2. 



APPENDIX D: BOUNDS ON E Q 

In order to derive the bound first consider the contribu- 
tion of the charged sector in the absence of the monopole 
antimonopole pair. The number of modes with frequency 
less than or equal to A + 1 is clearly 

N n J 

r, = EE E 2=-A^V + l)(2A + 7). (Dl) 

n=l J=l M=-J 

Let us denote the frequency of the T th mode by A(r). 
Then A(r) is is a staircase function that jumps at T = 
T N where A = 1,2,.... Evidently A(r) = N + 1 for 
r^r-i < r < r^v- In terms of this staircase function 
the contribution to E Q is given by (1/2) / o rc dTA(T). We 
can derive upper and lower bounds by approximating the 
staircase A(r) by smooth functions that intersect the top 
of each step and the bottom of each step respectively. 
Explicit calculation leads to the bounds 



RE U 



1 /3 



RE] 



4 \2 
1 /3 



4/3 



r y 3 + t Tc 



lower 



4 V 2 



4/3 



r 4 / 3 - -V 



(D2) 



(C7) These expressions apply in the relevant limit T c 3> 1. 
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The zero-point energy of the charged sector when the 
pair is present can be bounded similarly. We make the 
additional simplification that A n j = n whilst deriving 
the lower bound and A„j = n + 1 whilst deriving the 
upper bound. The results are: 



RE\ 

ower 



1 /3 
4 



1 /3 



4 V 2 



4/3 



r 4/3 _ i r 

4 c 



(D3) 



Upper and lower bounds on Eq may now be derived by 
taking the differences of the bounds in Eqs. (|D2jl and 
(|D3j) . In the discussion so far we have concentrated 
upon a single charged sector. Since there are in fact two 
charged sectors, + and — , we need to double the answer, 
leading to the result: 



T c < RE Q < 2T C 



(D4) 
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